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Abstract 

We develop a very simple method to study the high temperature, or equiv- 
alently high external field, behavior of the Sherrington-Kirkpatrick mean 
field spin glass model. The basic idea is to couple two difi'erent replicas 
with a quadratic term, trying to push out the two replica overlap from its 
replica symmetric value. In the case of zero external field, our results re- 
produce the well known validity of the annealed approximation, up to the 
known critical value for the temperature. In the case of nontrivial exter- 
nal field, we can prove the validity of the Sherrington-Kirkpatrick replica 
symmetric solution up to a line, which falls short of the Almeida-Thouless 
line, associated to the onset of the spontaneous replica symmetry break- 
ing, in the Parisi Ansatz. The main difference with the method, recently 
developed by Michel Talagrand, is that we employ a quadratic coupling, 
and not a linear one. The resulting flow equations, with respect to the pa- 
rameters of the model, turn out to be much simpler, and more tractable. 
As a straightforward application of cavity methods, we show also how to 
determine free energy and overlap fluctuations, in the region where replica 
symmetry has been shown to hold. It is a major open problem to give 
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a rigorous mathematical treatment of the transition to rephca symmetry 
breaking, necessarily present in the model. 



1 Introduction 

The mean field spin glass model, introduced by Sherrington and Kirkpatrick in 
[0, is here considered in the high temperature regime, or, equivalently, for a large 
external field. It is very well known, on physical grounds, that in this region the 
replica symmetric solution holds, as shown for example in 0, and references 
quoted there. However, due to the very large fluctuations present in the model, 
it is not so simple to give a complete, mathematically rigorous, characterization 
of this region, especially when there are external fields. Rigorous work on this 
subject include 0, ||5[, 0, [||. For other rigorous results on the structure 
of the model, we refer to g, ig, 0, H, IH, pf. 



The method developed in by Michel Talagrand, is particularly interesting. 
The starting point is given by the very sound physical idea that the sponta- 
neous replica symmetry breaking phenomenon can be understood by exploring 
the properties of the model, under the application of auxiliary interactions, which 
explicitly break the replica symmetry. In the replica symmetric solution is 
shown to hold in a region, which (probably) coincides with the region found in 
the theoretical physics literature, as shown for example in [0, i.e. up to the 
Almeida-Thouless line ||16 . 



The main tool in Talagrand's treatment is an additional minimal replica cou- 
pling, linear in the overlap between two replicas. Then, a kind of quadratic 
stability for the so modified free energy leads immediately, through a general- 
ization of the methods developed in to establish the validity of the replica 



symmetric solution in a suitable parameter region. 

Here we propose a very different strategy, by introducing a quadratic replica 
coupling, attempting to push the overlap away from its replica symmetric value. 
In a sense, our method is the natural extension, with applications, of the ideas 



put forward in |T^, where sum rules were introduced for the free energy, by ex- 
pressing its deviation from the replica symmetric solution in terms of appropriate 
quadratic fluctuations for the overlap. We choose exactly these quadratic fluctu- 
ation terms to act as additional interaction between two replicas, thus explicitly 
breaking replica symmetry. Then, a generalization of the sum rules, given in [^, 
for this modified model, allows us immediately to prove that the free energy of 
the original model converges, in the infinite volume limit, to its replica symmetric 
value, at least in a parameter region, explicitly determined. 

The organization of the paper is as follows. In Section 2, we recall the basic 
definitions of the mean field spin glass model, and introduce the overlap distri- 
bution structure. As a first introduction of our method of quadratic coupling, in 
Section 3, we treat the well known case of zero external field, by showing that the 
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annealed approximation holds, in the infinite volume limit, up to the true critical 
inverse temperature /3c = 1- Our proof shows explicitly that there is a strong 
connection between the critical value of the transition temperature for the zero 
external field model, and the analogous, and numerically equivalent, temperature 
for the well known ferromagnetic Curie- Weiss mean field model. 

In Section 4, we consider the model with external field, and introduce the as- 
sociated model with quadratic replica coupling. Then, simple stability estimates 
give immediately the convergence of the free energy to its replica symmetric value, 
in a suitable, well defined, region of the parameters. 

Section 5 reports about results on the free energy and overlap fiuctuations, in 
our determined replica symmetric region. We also sketch the method of proof, 
based on cavity considerations, as developed for example in and A more 



complete treatment will be found in a forthcoming paper [17 



Finally, Section 6 is dedicated to a short outlook about open problems and 
further developments. 

For the relevance of the mean field spin glass model for the understanding of 



the physical properties of realistic spin glasses, we refer to [18|, but see also [19 



2 The general structure of the mean field spin 
glass model 

The generic configuration of the mean field spin glass model is defined by Ising 
spin variables ai = ±1, attached to each site i = 1,2,...,N. The external 
quenched disorder is given by the N{N — l)/2 independent and identical dis- 
tributed random variables Jij, defined for each couple of sites. For the sake of 
simplicity, we assume each Jij to be a centered unit Gaussian with averages 

E(J,,)=0, Ei4) = l. 

The Hamiltonian of the model, in some external field of strength h, is given by 

Hn{(t, J) = --j^ JijCTiaj -h^^ai. (1) 

(id) « 

The first sum extends to all site couples, an the second to all sites. The nor- 
malizing factor 1/\/N is typical of the mean field character of the model, and 
guarantees a good thermodynamic limit for the free energy per spin, i.e., the 
existence of a finite and non trivial limit for the free energy as ^ oo. The 
first term in (|l|) is a long range random two body interaction, while the second 
represents the interaction of the spins with a fixed external magnetic field h. 

For a given inverse temperature f3, we introduce the disorder dependent parti- 
tion function Z]y{P, J), the (quenched average of the) free energy per site fN^P), 
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the internal energy per site UN{f3), the Boltzmann state cuj, and the auxihary 
function a^lP), according to the definitions 

Z^i(3,J)= J2 exp{-i3HNia,J)), (2) 

ITl...(TiV 

- (3fN{P) = N-'E\ogZ^{P, J) = aNiP), (3) 
ujiA) = Zjv(/3, J)-' Yl Aexp{-f3HNia, J)), (4) 

(T1...(T]V 

un{(3) = N-'Eu;j{HN{a, J)) = dp{(5fNm = -a^«iv(/9), (5) 

where A is a generic function of the cr's. In the notation cuj, wc have stressed the 
dependence of the Bohzmann state on the external noise J, but, of course, there 
is also a dependence on /?, h and A^. 

We are interested in the thermodynamic limit N ^ oo. 

Let us now introduce the important concept of replicas. Consider a generic 
number s of independent copies of the system, characterized by the Boltzmann 
variables al^\ af^, . . ., distributed according to the product state 

O -,,(1), ,(2) 

where all cJj"^ act on each one (T^^^^'s, and are subject to the same sample J of the 
external noise. Clearly, the Boltzmannfaktor for the replicated system is given 

by 

exp {-P{Hr,{a^'\ J) + Hr,{a^'\, J) + . . . + Hi,{a^'\ J))) . (6) 
The overlaps between two replicas a, b are defined according to 

i 

and they satisfy the obvious bounds 

-1 < Qab < 1. 

For a generic smooth function F of the overlaps, we define the () averages 
{F{qu,qn,...)) = EnjiF{qu,qn, ■■■)), 

where the Boltzmann averages Qj acts on the replicated a variables, and E is 
the average with respect to the external noise J. 

We remark here that the noise average E introduces correlations between 
different groups of replicas, which would be otherwise independent under the 
Boltzmann averages as for example qi2 and ^34. 

The averages are obviously invariant under permutations of the replicas. 
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Overlap distributions play a very important role in the theory. For example, 



a simple direct calculation shows [g , |TT[ that 



5/3«iv(/5) = f(l-(gL))- (7) 

In order to introduce our treatment, based on flow equations with respect to 
the parameters of the theory, it is convenient to start from a Boltzmannfaktor 
given by 

exp 1 ^ JijCTiaj + (3hY^ai + y/x^^JiaA . (8) 

Notice that we have introduced an auxiliary additional one body random in- 
teraction ruled by the strength y/x, x > 0, and quenched independent and 
identically distributed centered unit Gaussian random variables Jj, so that 

In order to get the original model, we have to put a; = at the end. Moreover, 
we have written P = y/t, with t > 0. The variables t, and x, will play the role 
of time variable, and space variable, respectively, in our flow equations. Now 
we deflne the partition function Z by using the Boltzmannfaktor (H), and the 
auxiliary function ajv(x,t) in the form 

aN{x,t) = N-^ElogZ^. 

Then, as in the proof of (^, we have 

dtar,{x,t) = ^il-{ql^)) (9) 

d,aNix,t) = ^{l-{qi2)). (10) 

It is very simple to calculate explicitly the average N'^Elog Zn for t = 0, at 
a generic strength xq of the one body random interaction. In fact, at t = 0, 
the interaction factorizes, and the spins at different sites become independent. 
Therefore we have 



a(xo, 0) = log2 + J logcosh.{Ph + z\^Q)dfj.{z), 



(11) 



independently of A^, where dfi is the centered unit Gaussian, representing each 
single Ji, 



djj{z) = exp{——)dz/\'2TT. 
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Starting from ([Tl|), (|T0|) at t = 0, we can immediately calculate the order param- 
eter g(xo) according to 

Qi^o) = (^12) (3^0,0) = J tanli^ {Ph + Zy/x^)dfi{z). 

Let us now consider consider linear trajectories at constant velocity given by 

x{t) = Xq- q{xo)t, (12) 

where Xq is the initial starting point. It is immediate to verify that there is only 
one trajectory passing for a generic point {x,t), x > 0,t > 0, i.e. the initial point 
is uniquely determined by {x,t), and so is for q. In fact, the following Theorem 
holds. 

Theorem 1. Consider first the case of non zero external field h. Then, for a 
generic point {x,t), x > 0, t > 0, there exists a unique Xo{x,t) such that 

X = Xoix,t) -q{xo{x,t))t, 

and a unique q{x,t) = q{xQ{x,t)), such that 

q{x,t) = / tanh^(/5/i + z^y x + q{x, t)t)dfi{z). 



Ifh = 0, then invertibility is assured, with some XQ{x,t) > 0, in the region x > 0, 
if: > 0, with the exclusion of the segment x = 0, < t < 1. 

The proof is very simple, and can be found in . 



By following the methods of |]T5|], sum rules connecting a and its SK approx- 



imation can be easily found by using transport equations. They involve overlap 
fluctuations. In fact, let us consider a along the trajectories in (|T2|), i.e. a{x{t),t). 
An easy calculation gives 

f^a^{x{t),t) = \{l-q{xo)f - \{{qu-qr). (13) 

Now we integrate along t, take into account the initial condition (|11|), and the 
invertibility of (0), and find the sum rule 

1 /■* _ 

a{x,t) = aNix,t) + - {{qi2~ qY)x{t'),t'dt' . (14) 
4 Jo 

Here, we have defined the replica symmetric Sherrington-Kirkpatrick solution |1|], 
0, in the form 



a 



(x, t)=log2 + j \ogcosh.{l3h + z^Q)d^{z) + - q{xQ)f . (15) 



A very simple, but important, consequence of the sum rule is that is domi- 
nated by its replica symmetric solution, uniformly in N , 

aNix,t) = N'^ElogZ^ <a{x,t). (16) 

This is a simple consequence of the posit ivity of the term under integration in 

Now, we are ready to explain our method of quadratic coupling, starting with 
the simple case of zero external field, and then going to the case of nontrivial 
external fields. 



3 Quadratic coupling for zero external field 

The high temperature region (/3 < 1) of the zero external field SK model is a 
very particular case where everything can be computed. As it is well known 0, 
in this case the annealed approximation is exact in the infinite volume limit. In 

fact, we have 

jim l^ln Z^{t, J) = ait) ^ ln2 + ^ = ^ lnEZ^(t, J) + ^. (17) 

In this section we give a new proof of Eq. (|I7p, based on sum rules for the 
free energy. Our method is very simple and can be easily extended to the case of 
nontrivial external field, considered in the next Section. 

When h = 0, X = 0, then also q = and Xq = 0. Then, the sum rule in (0) 
reads ^ 

aN{t) = a{t)-^j^{ql,)t>dt'. (18) 

The presence of (^12)' order parameter, suggests to couple two replicas with 
a term proportional to the square of the overlap, the corresponding partition 
function being 




with A > 0. The effect of the added term is to give a larger weight to the 
configurations having qi2 7^ 0, thus favoring nonselfaveraging of the overlap. Of 
course, the system possesses spin-flip symmetry also for A 7^ 0, so that (^12) = 
0. Therefore, if (g'12) 7^ then the overlap is nonselfaveraging. Now replica 
symmetry is explicitly broken. 

The basic idea of our method is to show that, as long as t < 1 and A is 
small enough, the term \Nq\2 does not change the value of the free energy in the 
thermodynamic limit. Therefore, "most" configuration must have qi2 = and 
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the overlap must be self aver aging. In order to implement this intuitive idea, one 
introduces the A dependent auxiliary function 

aN{t,\) = In Zn, 

where the normalization factor 1/2 is chosen so that a]sf{t,0) = aAr(t). Through 
a simple explicit calculation, we can easily calculate the t derivative in the form 

dtaN = ^ (1 + (??2) - 2(g?3)) , (20) 

where now all averages ( ) involve the A-dependent state with Boltzmannfaktor 
given in agreement with (|TPp. Moreover, it is obvious that 

dx&N = ^{qu)- 

Starting from some point Aq > 0, consider the linear trajectory A(t) = Xq — t, 
with obvious invertibility in the form Aq = A + t. Let us take the t derivative of 
cun along this trajectory 

^a^(t, A(t)) = {dt -dx)aN = ^- ^(g?3)t,A(t)- 

Notice that the term containing (g^g) disappeared. By integration we get the sum 
rule and the inequality 

t ^ 1 /"* t ^ 

ttNit, A) = - + aAr(0, -^o) - 2 J (^13)t',A(i') dt' < - + aAr(0, Ao), (21) 

where {q\^)t' ,\{t') refers to \{t') = Xq — t' = X + t — t' . 

Next, we compute (5Ar(0, Aq). We introduce an auxiliary unit Gaussian z, and 
perform simple rescaling, in order to obtain 

«iv(0,Ao) = ^lnE^^'°''''^ = ^l^E I e^''^' d^,{z) (22) 

{cr,o-'} {a,o-'} 

= ^^2 + ^lny" fcosh j d^i{z) 

= ln2 + ^ln I rfyy^expiv('-Ao^ + lncosh(yAo)] ,(23) 



where we performed the change of variables z to yy/NXo in the last step. It is 
immediately recognized that the integral in (^) appears in the ordinary treat- 
ment of the well known ferromagnetic mean field Curie- Weiss model. The saddle 
point method gives immediately 

1 / y2 \ 

lim a;Ar(0, Aq) = ln2 H — max — Aq h lncosh(?/Ao) . (24) 

N^oo 2 y \ 2 / 
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Therefore, the critical value for Aq is Ac = 1. For Aq > 1 we have 

hm aAr(0, Aq) > In 2, 

while for Aq < 1, one can use the elementary property 2 In cosh x < x"^ to find 

«^(0,Ao)<ln2 + ^ln^^. (25) 

Notice that, when Aq approaches the value 1~, the term of order diverges, 
since Gaussian fluctuations around the saddle point become larger and larger. 
Thanks to (25), the inequality in ( pi]) becomes 

ajv(t, A) < a(t) + -^In- 



4N 1 - Ao 

which holds for < Aq < 1, i.e., for < t + A < 1. 

Next, we use convexity of aN{t, A) with respect to A and the fact that 

to write ^ ^ ^ 

(^N{t) + -{ql2)t <«7v(t,A) < + 

for A > 0. For < t < t < 1, choose A = (1 - 1)/2, so that 

A + t<Ao = (l + t)/2<l, 

and 



< \iait) - a^it)) + ^ln_. (26) 

Recalling Eq. ([I8|), one has 

d XI XX 

-(«(t) - aM{t)) = -{ql,)t < j{a{t) ' a^W) + Jf;^^^Y^o' ^^^^ 

so that 

aN{t) = a{t) + 0{l/N), 

uniformly for 0<t<t<l. Of course, from Eq. (PT| ) and convexity of one 
also has 

a^(t,A) =a{t) + 0{l/N), 
{Qlskx = 0{1/N), 

for < t + A < Ao < 1. 

We have gained a complete control of the system in the triangular region 
0<t<l,0<A<l — t. Note that we have not only proved Eq. (|17D but we 
also shown that the leading correction to annealing is of order at most 
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4 The general case 

The method we follow for the general case, where the Boltzmannfaktor is given 
by (j^), is a direct generalization of the one explained in the previous section. In 
fact, by taking into account the t derivative in (p!3D, we are led to introduce the 
auxiliary function 

aAr(x, A,t) = — ElnZiv(x, A,t; J), 

where Zn is the partition function for a system of two replicas coupled by the 
term 

with A > 0. In order to simplify notation, we omit the argument h. 
Now the t derivative is given by 

1 A 

dtUN = 4(1 + (912) - 2(gi3)) + 2 (9 - (512))^, 

while the x and A derivatives appear as 

1 A dq 

d^&N = 2 + ^^12) - 2(gi3)) + -{q - (^12))^, 

dx&N = ^ {{qi2 - q{x,t) f) ■ 

Starting from points Aq > 0, xq, consider the linear trajectories X(t) = Xq — t, 
x{t) = xo — q{xo)t, as in (^2|), with obvious invertibility as explained before. Since 
q is constant along the trajectory, one finds for the total time derivative of 

j^aN{x{t),\{t),t) = {dt- qdx - dx)dN = ^(1 - g)^ - ^ ((gis - g)^> • 

Notice that in this case the term containing {{qi2 — q)^) disappeared. 
By integration, we get the sum rule and the inequality 

a^(x,A,t) = -(l-g)^ + aiv(a;o,Ao,0)-- / ((fe - ^)^>t,^^(4,)^A(t') 

^ 

< ^(l-g)^ + a7v(xo,Ao,0), 

where ((gis - q)^)t>^^(^t')Xt') ^^^^^^ *° 

A(t') = X^-t' = X + t-t', 
x{t') = Xq — qt' = X + q{t — t'). 
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If f2j is the product state for two replicas with the original Boltzmannfaktor given 
by (H), then we can write 

aN{x, A, t) - aAr(x, t) = —E In Vtj \ exp -AiV(gi2 - qf 

Therefore, by exploiting Jensen inequality, we have, for A > 0, 

^ ((912 - qf)^^t < aN{x, A, t) - aNix, t). 
Let us also define 

An{xo, Ao) = a7v(xo, Aq, 0) - a{xo, 0) = In n^j ^exp ^AoA^(gi2 - g)^^ , 

(28) 

where we have introduced the state Q^j for two replicas, corresponding to t = 0, 
and X = Xq, in (j^). Notice that f2j is a factor state over the sites i. 

By collecting all our definitions and inequalities, and taking into account the 
definition (plSf ), we have 

^ ((^12 - Qy)^^t < Aiv(xo, Ao) + a{x, t) - aN{x, t). 
Let us now introduce \c{xq) such that, for any Aq < Ac(xo), one has 

lim A7v(xo, Ao) = 0. 

Then, by the same reasoning already exploited starting from (pG]), and taking 
into account (p!3D , we obtain the proof of the following 

Theorem 2. For any t < Ac(xo(x,t)), where xo{x,t) is defined as in Theorem 
d, we have the convergence 

lim aN{x,t) = d{x,t). (29) 

N—*oo 

For the specification of Ac(xo), we can easily establish the complete charac- 
terization of the Ajsf limit. In fact, the following holds. 

Theorem 3. The infinite volume limit of An is given by 

lim AAr(xo,Ao) = A(xo,Ao). 

Here, A(xo, Ao) is defined through the variational expression 

A(a;o, Ao) = - max ( / ln(cosh fi + tanh^(/3/;. + z^/xo) sinh fi)dp{z) — fiq — -^r- 
2 tJ. \J 2Ao 

(30) 

where dp{z) is the centered unit Gaussian measure. 
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Of course, the expression ( pOD is in agreement with (0), when there are no 
external fields. 

It is easy to realize that the value Ac, in the general case, is strictly less than 
the expected value tc, following from the Almeida-Thouless argument. 



tc J cosh '^(Ph + z^/xo)dfl{z) = 1. 



The proof of the Theorem ^ is easy. First of all let us establish the elementary 
bound, uniform in N, 

A7v(xo,Ao) > A(xo,Ao). (31) 

In fact, starting from the definition of Ai\f{xo, Aq) given in (^), we can write, for 
Ao 7^ 0, and any /i, 

\ 2 



(gi2 -g) > 2— (gi2 -g) - . , 
and conclude that 

AAr(xo, Ao) > aoifi) - — , (32) 

where we have defined 

ao(/i) = In (exp ^iV(gi2 - q)) 

If 1 

= - / In (cosh /i + tanh^ {(3h + z^/xo) sinh /x) dp{z) fiq. 

2, J 2 



Of course, it is convenient to take the max^ in the r.h.s. of (|32D , so that the bound 
in (|3T|) is established. The proof that the bound is in effect the limit, a.s N —>■ oo, 
can be obtained in a very simple way by using a Gaussian transformation on (pSf) , 
as it was done in (|22|). In fact, we now have 

In (^exp hoNiq,, - q)'^ (33) 

= ^E\n J fiO (^exp yVV(gi2 - q)z) dp{z) 

Therefore, by exploiting the fact that also fij factorizes with respect to the sites 
i, we can write 

A^(xo,Ao) = -^E\n j JJ {^osh ^z + iaTih\(3h + Ji 

exp {-^JXoNqz^ dp{z). (34) 



xo) sinh \l ^z 




X 
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Now, we find convenient to introduce a small e > 0, so that 



11 

3- = T7 + e. (35) 

Notice that Aq < Ag. We also introduce the auxihary (random) function 

1 1 
<pNiy, Aq) = X]^^ (coshy + tanh^(/?/i + Ji^/x^) sinhy) -^V- 

i " 

By the strong law of large numbers, as — > 00, for any we have the J almost 
sure convergence of (pN{y, Aq) to (f){y, Aq) defined by 



0(?/, Aq) = J \n (^cosh y + tanh^{P h + Zy/lEo) sinli dp{z) — qy — 
= E(f)N{y,X'f^). 



Ill 
2A'o 



Here dp{z) performs the averages with respect to the Ji variables. Let us also 
remark that the convergence is J almost surely uniform for any finite number of 
values of the variable y. 

Now we can go back to (0), write explicitly the unit Gaussian measure dp{z), 
perform the change of variables y = zy/\oN~^, make the transformation (|35D , 
take the sup^^ for the (j)N, and perform the residual Gaussian integration over y. 
We end up with the estimate 

An{xo, Ao) < sup (pN^V, Aq) + 7^ In -^=. (36) 
2 y 2I\ VAoe 

Since the J-dependent sup^ is reached in some finite interval, for any fixed Aq, 
and the function 0Ar is continuous with respect to y, with bounded derivatives, 
we can perform the sup^ with y running over a finite discrete mesh of values, by 
tolerating a small error, which becomes smaller and smaller as the mesh interval 
is made smaller. But in this case the strong law of large numbers allows us to 
substitute (pN with 0, in the infinite volume limit — > 00. On the other hand, 
the second term in the r.h.s. of (^) vanishes in the limit. Therefore, we conclude 
that ^ 

limsup AAr(xo, Ao) < - sup </>(?/, Aq). 

Af— >00 2 y 

From continuity with respect to Aq, we can let Aq approach Aq, and the Theorem 
is proven. 



5 Fluctuations of overlaps and free energy 

In the previous sections we proved that, in a certain region of the parameters 
t,x,(3h, the typical values of the free energy In Z^/N and of the overlap qab are 
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the replica symmetric expressions a and g, respectively. In the same region, one 
can obtain a more precise characterization of the fluctuations of these quantities, 
for iV — cxo, showing that a central limit-type theorem holds, after suitable 
rescaling. This will be analyzed in detail in a subsequent paper |T^. Here, we 
just give the main results and sketch the ideas underlying the proof. 

Concerning the fluctuations of the overlap around the Sherrington-Kirkpatrick 
order parameter q, we prove the following 

Theorem 4. j7^/ The rescaled overlap variables 

d = ^{Qat - q) 

tend in distribution, as N oo, to centered, jointly Gaussian variables with 
covariances 

{ej = A{t,x,h) 

(LbLc) = B{t,x, h) 

i^ab^cd) = C(t,X, h), 

where b ^ c, c a,b and d ^ a,b. The expressions of A, B, C are explicitly given 
and coincide with those found in /j7^. 



Recently, an analogous result was proved independently by Talagrand [ pO|] , 
who computed the N —>■ oo limit for all moments of the ^ variables. 

The scheme of our proof is as follows: The control we obtained on the coupled 
two replica system and concentration of measure inequalities for the free energy 



2^ imply that the fluctuations of qab from q are exponentially suppressed for N 
large. Then, by means of the cavity method [@] one can write a self-consistent 
closed equation for the characteristic function of the variables This equation 
turns out to be linear, apart from error terms which vanish asymptotically for 

^ oo, thanks to the strong suppression of the overlap fluctuations. The 
solution, which is easily found, coincides with the characteristic function of a 
Gaussian distribution, with the correct covariance structure. 

Concerning the free energy, the result we prove is the following: 



Theorem 5. / [7^/ Define the rescaled free energy fluctuation as 



Then, 



r f. I, n- /.r i'^'^ZN{t,x,h;J) _ 
fN{t, x,h;J) = VN[ — a{t, x, h) 



fN(t, X, h] J) Af{0, cr^(t, X, h)), 



where N'{m, a^) denotes the Gaussian random variable of mean m and variance 
cr^; and 

a'^it, X, h) = Var (in cosh.{zy/qt + x + ph)^ —. 

Here, Var{.) denotes the variance of a random variable and z = A/'(0, 1). 
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This result is a consequence of Theorem ^ and of concentration of measure 
inequahties for the free energy. 



6 Outlook and conclusions 

We obtained control on the thermodynamic limit of the model, in a region above 
the Almeida- Thouless line, by suitably couphng two replicas of the system and 
studying stability with respect to the coupling parameter. The question naturally 
arises, whether and how this method can be extended up to the expected critical 
line. This problem seems to be common to all approaches proposed so far. 

The method can be also further generalized to the case where more and more 
replicas are mutually coupled. In this case, replica symmetry is explicitly broken 
at various levels, and it is possible to give a generalization of the Ghirlanda- 



Guerra relations [13|. We plan to report soon on these generalizations |22 
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